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CONVERSION FACTORS, U.S. CUSTOMARY TO METRIC (SI) 
UNITS OF MEASUREMENT 


U.S. customary units of measurement used in this report can be converted to metric (SI) 


units as follows: 


Multiply by To obtain 
inches 25.4 millimeters 
2.54 centimeters 
square inches 6.452 square centimeters 
cubic inches 16.39 cubic centimeters 
feet 30.48 centimeters 
0.3048 meters 
square feet 0.0929 square meters 
cubic feet 0.0283 cubic meters 
yards 0.9144 meters 
square yards 0.836 square meters 
cubic yards 0.7646 cubic meters 
miles 1.6093 kilometers 
square miles 259.0 hectares 
acres 0.4047 hectares 
foot-pounds 1.3558 newton meters 
ounces 28.35 grams 
pounds 453.6 grams 
0.4536 kilograms 
ton, long 1.0160 metric tons 
ton, short 0.9072 metric tons 
degrees (angle) 0.1745 radians 
Fahrenheit degrees 5/9 Celsius degrees or Kelvins! 


' To obtain Celsius (C) temperature readings from Fahrenheit (F) readings, use formula: C = (5/9) (F — 32). 
To obtain Kelvin (K) readings, use forumla: K = (5/9) (F — 32) + 273.15. 


STORM SURGE SIMULATION IN TRANSFORMED COORDINATES 


Volume II. Program Documentation 


by 
John J. Wanstrath 


I. INTRODUCTION 


Presented in this volume are the program documentation and list- 
ings of the coded programs for carrying out a simulation of a 
hurricane-induced storm surge on the Continental Shelf in curvilinear 
coordinates for a given reach of shelf. There are three separate 
programs detailed here for this purpose. 


The first program, CONFORM, is employed for evaluation of the 
transformation coefficients which map the given reach of the Conti- 
nental Shelf into a rectangle in the image plane, the shelf region 
being delineated by a smoothed version of the coastline, the shelf 
break (taken as the 180-meter depth contour in the examples), and 
bounded laterally by two parallel straight lines. The theory and 
several example applications of the transformation procedure are 
given in Section II of Volume I. The example input for CONFORM con- 
tained in the listings here are for the mapping of the region from 
a section across Laguna Madre about 90 kilometers south of Brownsville 
to Marsh Island. Particular care must be taken to follow the pro- 
cedure exactly for the example if CONFORM is to be verified by obtain- 
ing the given transformation coefficients. This procedure is given 
explicitly in the CONFORM documentation. 


The second program, GRID, develops the detailed computing grid 
information, based on the mapping coefficients evaluated by CONFORM 
plus coordinate stretching information supplied by the user (see 
Section III of Volume I). Part of the output of GRID is a listing of 
the grid positions which are required by the user in order to read 
from appropriate charts the detailed bathymetry field, which is 
necessary input for the final program SSURGE. The example data 
supplied here for GRID are for the Laguna Madre to Marsh Island 
region. 


Program SSURGE (Shelf Surge) carries out the numerical integra- 
tion of the storm surge equations in the transformed coordinate system 
supplied by CONFORM and GRID, using a parametric representation of a 
hurricane wind field and pressure field. The theory is given in Sec- 
tion III of Volume I. The particular example data given here are for 
Hurricane Carla and the Laguna Madre-Marsh Island grid system. 


The appendixes to this volume contain detailed FORTRAN listings 
of the three programs in this application. The data to be supplied 
by the user for other applications are discussed in the documentation 
of each of these programs. 


II. COMPUTER PROGRAM DOCUMENTATION FOR PROGRAM-CONFORM 


1. Program Purpose. 


The purpose of the program is to determine the transformation 
coefficients which will conformally map the interior region bounded 
by the actual coastline, a seaward boundary curve, and two parallel 
lateral boundaries into a rectangle in the image plane. 


2. Program Description. 


The program is written in FORTRAN IV language. This program and 
the program GRID provide all the necessary computing grid data for 
input to program SSURGE. The program GRID takes, as input, the 
transformation coefficients and determines the computing grid infor- 
mation (such as, scale factors, grid point locations, and, at each 
grid point, the orientation of the &-axis to the x-axis). 


The program CONFORM is composed of: 


MAIN Defines constants. Reads and writes the 
coordinates delineating the given coastline 
and seaward boundary curve. Calls Subroutine 
COEFFS. 


SUBROUTINE COEFFS Determines the transformation coefficients. 
At the completion of each iteration, the 
coefficients, the variance between the 
transform-generated curves and that Speci- 
fied, and other pertinent information are 
written. COEFFS interfaces all other 
program subroutines and functions. 


FUNCTION XTRAN Is the transformation function x(&,n) 
FUNCTION YTRAN ' Is the transformation function y(é,n) 
SUBROUTINE SLFAC Determines the scale factor and derivatives, 


ox/9—& , dy/d& , for a given value of €& and 
m6 


SUBROUTINE CURLYB Determines the necessary parameters to fit a 
spline under tension through the given coast- 
line coordinates. The spline is fitted with 
Y2 as a function of X2 


FUNCTION CURVYB Interpolates the given coastline, returning 
a value of y at a specified value for x. 


FUNCTION CURDYB Differentiates the given coastline, returning 
a value of dy/dx at a specified value for 
xe 


SUBROUTINE CUR2YB Determines the necessary parameters to fit a 
spline under tension with X2 as a function 
of coastline arclength. 


FUNCTION CUR4YB Interpolates the given coastline returning a 
value of x at a specified value for arclength. 


SUBROUTINE CURSYB Determines the necessary parameters to fit a 
spline under tension with Y2 as a function of 
coastline arclength. 


FUNCTION CURSYB Interpolates the given coastline returning a 
value of y at a specified value for 
arclength. 


There are identical subroutines and functions as delineated above 
for the seaward boundary curve specified by coordinates X2P and 
Y2P . These subroutines and functions are recognized by the same 
names as their counterparts with a terminal letter A. For example, 
SUBROUTINE CUR1YA determines the necessary parameters to fit a spline 
under tension through the given seaward boundary curve. 


3. Type of Computer, 


The program CONFORM can be run on any computer with minimum core 
requirements of approximately 24K (based on the present sample program). 
However, significantly more computer memory would be required if one 
desires a large number of coefficients and/or numerous integration 
points, The program has been executed successfully on IBM 360, CDC/ 
6600 and 7600, and GE/635, The present sample program requires no 
auxiliary storage devices, peripherial devices, or magnetic tape input 
or output. No site-orientated computer plot. routines are involved in 
the program, Approximately 20 minutes of machine time on a CDD/7600 
is required for the sample program (total number of coefficients, 

2 x NMAX = 220; number of integration points, 0 < —€ <r , = 110; and 
number of iterations = 80). 


4, Input Data, 


Input data, other than constants defined in MAIN, are read in 
MAIN on IBM cards prepared according to the following list: 


Gl) Gases i 
IWANT, MQ, MOP, NMAX, Continuation flag, number of shore- 
JMAS1, IL, VARWT line and seaward boundary coordinates 


number of mapping coefficients, maximum number of 
iterations, number of integration points, and the 
convergence criterion in format 615, F5.3. 


(2) Card Group 2 


NZ NZ The x,y coordinates (units in x,y space) of 
the given coastline in the region 0 <x<A 
are read with one pair per card in format 3X, 
B76 Sho 764 (laine i'50)) . 


(3) Card Group 3 


X2P,Y2P The x,y coordinates (units in x,y space) of 
the given seaward boundary curve in the region 
0 <x < are read with one pair per card in 
format F7.2, 3X, F7.2 (limit 150). 


Optional Card 4, Card Group 3, and Card 6 


If IWANT = 1, indicating the program is being re- 
started, the following cards must be supplied: 


(4) Card 4 


B, BZRO whe values of s)8ayand s Bossi units of length of 
x,y space from the last iteration of the prev- 
ious run in format 2E14.7. 


(S) Card Group 5 


COB, COC The NMAX cards containing the dimensionless 
Fourier-type transformation coefficients from the 
last iteration of the previous run in sequential 
order with one pair per card (format 2E14.7). If 
more coefficients are desired in the present run 
than the previous one, blank cards should be 
supplied for the difference. 


(6) Card 6 


WSX, WSY The value of the weighting factors for the sea- 

WCX , WCY ward boundary curve (x and y component) and 
coastline (x and y component) from the last 
iteration of the previous run (4E14.7). 


5. Mathematical Procedures and Program Limitations. 


The conformal mapping relations, solutions for the transformation 
coefficients, and the iterative procedure for determining the coeffi- 
cients are presented in Volume I of this report. The mapping equa- 
tions are sufficiently general to treat the situation where either or 


both given curves are multivalued in y for a specified x value. 
The only program limitation relates to the computer memory capacity. 


No program error messages or consistency checks are incorporated in 
this routine. 


Ge RlownGhantr 


MAIN - CONFORM 


Read constants; MQ, MQP, NMAX, 
IL, JMAX1, VARWT, IWANT 


Read and write coordinates delineating 
the coastline and seaward boundary curve; 
X2 WA goal XA 5 AAP 


a 


inaktevauleezen Ay bee ROMs m NS X 
WCX, WCY , COB(N) and COC(N) 


Read transformation 
coefficients and 
weighting factors 


CALL SUBROUTINE COEFFS 


SUBROUTINE COEFFS 


Compute coastline arclength 


Arclength computed directly 
from given coordinates, 
O25 2) 

Yes Consequently, assumed that 
used has sufficient data 
points to accurately define 
the curve. 

Stored in ARCYB 


COASTLINE HAS 
FOLDBACK, 
TYB=1 


No 


_ CALL CUR1YB 
Arclength computed from a spline 
fitted through Y2 as a function 
of X2 with 9 additional points 
between the given coordinates to 
further resolve the curve. 
Stored in ARCYB. 


CALL CUR2YB 
Fit a spline through X2 as a function of ARCYB 


CALL CUR3YB 
Fit a spline through Y2 as a function of ARCBY 


Compute seaward boundary arclength; the same 


procedure as for the coastline is followed. 


Arclength computed 
directly from X2P, 
WAP. 

Stored in ARCYA 


CALL CURLYA 
Arclength computed from a 
spline. Stored in ARCYA 


CALL CUR2YA 


CALL CUR3YA 


Write coastline and seaward boundary 
coordinates and arclengths. 


Compute CS and SN _ where 


CS = cos (nké) , 
: for n=1,2 --+ NMAX 
SN = sin (nk ‘ Zl 
ee ema ES(iet) Gite) 
where 
dai sAi oo Wb/A.)\ MWh ChyGine 
or 
i=1,2 °° IL/2+1, IL odd. 


Determine the x and y coordinates of the given 
coastline and seaward boundary curve from the arc- 
length of the transform-generated curves. 

ATXB , ATBY and ATXA , ATYA 


Compute B and _ BZRO 


Compute COB and COC 
for n=1,2 -- NMAX 


Write transformation coeffi- 
cients, iteration number, 
NMAX and IL. 


Compute the error function and variance between the trans- 
form-generated curves and that specified. Write the values. 
The x and y coordinates of the given curves are deter- 
mined from the arclength of the new transform-generated 
curves for the next iteration. Write the values. 


COMPUTED VARIANCE 


ICOUNT = ICOUNT + 1 


7. Glossary and Description of Terms. 


Arrays 

X2 Dimensioned MQ. Contains values of the x coordinate in 
units of x,y space of the given coastline (max, 150). 

WZ Dimensioned MQ. Contains values of the y coordinate in 
units of x,y space of the given coastline (max. 150). 

X2P Dimensioned MQP. Contains values of the x coordinate in 
units of x,y space of the given seaward boundary curve 
(max, 150). 

WiZB. Dimensioned MQP. Contains values of the y coordinate in 
units of x,y space of the given seaward boundary curve 
(max. 150). 

COB Dimensioned NMAX. Contains values of the dimensionless 
transformation coefficient BY (max. 200). 

COC Dimensioned NMAX. Contains values of the dimensionless 
transformation coefficient ce (max. 200). 

DUMB Dimensioned NMAX. Contains the iterative value for BO 
(max, 200). 

DUMC Dimensioned NMAX. Contains the iterative value for Co 
(max. 200). 


ARCYB Dimensioned MQ. Contains values of the arclength for the 
given coastline coordinates X2 , Y2 (max.150). 


ARCYA Dimensioned MQP. Contains values of the arclength for the 
given seaward boundary coordinates X2P , Y2P (max.150). 


SN Dimensioned NMAX x IL1 where IL1 is IL/2 for even IL or 
IL1 is IL/2 + 1 for odd IL. Contains values of sin (nké). 

CS Dimensioned NMAX x IL1. Contains values of cos (nké). 

ATXB Dimensioned IL. Contains values of the x coordinate of 


the given coastline as determined from the arclength of 
the transform-generated curve (max. 400). 


ATYB Dimensioned IL. Contains values of the y coordinate of 
the given coastline as determined from the arclength of the 
transform-generated curve (max. 400). 


ATXA 


ATYA 


Dimensioned IL. Contains values of the x coordinate of the 
given seaward boundary curve as determined from the arclength 
of the transform-generated curve (max, 400). 


Dimensioned IL. Contains values of the y coordinate of 
the given seaward boundary curve as determined from the 
arclength of the transform-generated curve (max.400). 


Dimensioned IL. Contains values of & (max.400). 


? Dimensioned IL. Temporary storage. 


Constants 


MQ 


MQP 


NMAX 


IL 


IWANT 


B 


BZRO 


WSX, WSY 


WCX, WCY 


Number of coordinates delineating the given coastline 
(max.150). 


Number of cgordinates delineating the given seaward boundary 
curve (max.150). 


Number of transformation coefficients, BO or Cc. (max. 200). 
Number of equally spaced integration points for 0 < € < xX 
with IL > NMAX (max. 400). 


If IWANT = 1, the program is to be re-started requiring 
input from the previous run. If IWANT # 1, it is the 
initial mapping of the region. 


Desired variance (in units of x,y space squared) between 
the transform-generated curves and that specified. . 


h = X2(MQ) = X2P(MQP) 
k = a/d 

-B 

B 

BS 


The x and y component of the weighting factor for the 
seaward boundary curve. 


The x and y component of the weighting factor for the 
coastline. 


8. Input and Output. 


The input data required by CONFORM to determine the transforma- 
tion coefficients which conformally map the Laguna Madre, Mexico, to 
Marsh Island, Louisiana, region into a rectangle are presented here 
as an example, The first card image gives the continuation code, 
number of coastline and seaward boundary coordinates, number of coef- 
ficients desired, maximum number of iterations, number of integration 
points, and the convergence criterion, The next 47 paired numbers 
are the coastline coordinates and the last 40 card images give the 
seaward boundary coordinates, 


Program CONFORM provides detailed and voluminous output concern- 
ing primarily the rate of convergence. These output statements are 
not necessary for program completion and can be easily deleted with 
little alteration to the sequence of instructions. Optional output 
statements are indicated in the program listing (App. A) by an arrow 
(<) on the right-hand side of the page. 


The results from CONFORM required by Program GRID are 8, B_ , 
and the coefficients, B, and C, . These are given in the next 
section as input to program GRID. 


Reference to Table 4 in Volume I of this report shows the values 
for the number of coefficients, maximum number of iterations, and 
number of integration points are only indicative of the final steps 
of this particular application. Explicitly, in order to obtain the 
given coefficients to conformally map the Laguna Madre to Marsh Island 
region into a rectangle in the image plane, the following steps must 
be followed: 


(1) Set NMAX to 40 and IL to 80 for the first 20 iterations, 
(2) For the next 10 iterations, NMAX is 60 and IL is 120, 


(3) Set NMAX to 80 and IL to 160 for iterations 31 through 
50. 


(4) Take NMAX as 90 and IL as 180 for the next 10 
iterations, 


(5) For iterations 61 through 70, NMAX is 100 and IL is 200. 
(6) Over the next 10 iterations, NMAX is 110 and IL is 220. 
(7) Starting with the 8lst iteration, WCX amd WSX are set 


to 0 and NMAX is 110 and IL is 220 through the 100th 
interation. 


(8) From iteration 101 through 110, NMAX is 130 and IL is 
260, 


(9) Over the last 10 iterations, NMAX is 150 and IL is 300. 


Note that steps (7) through (9) utilize the alternate solution 
to the mapping equations (9) and (10) in Section II of Volume I of 
this report, The instructions which must be altered or removed are 
indicated in the program listing by parenthesis with the proper 
instruction enclosed. 


The following input is required by CONFORM to conformally map the Laguna Madre, Mexico 
to Marsh Island, Louisiana region into a rectangle in the image plane. Note the fourth and 
fifth parameters on the first card image are only indicative of the final results. See 
Input and Output text for explanation of procedure used to obtain the mapping coefficients 
given as input to Program GRID. 


Cc 

Cc INPUT OATA 

Cc 
1 47 40 150 1500¢001 360200 110.00 

(2 c 

Cc Cc 

Cc COASTLINE COORDINATES Cc SEAWARC ECUNDARY COORDINATES 

Cc Cc 

Cc Cc 
000e00 C35e00 000-200 0CzZe00 
005e60 040200 007¢50 OCE&e0C 
010.50 044.80 016200 010200 
018.200 C5100 021-250 014-50 
028200 CSSe00 028200 0190200 
040.200 O0€Ee80 036200 022200 
050e50 C730e00 045-200 028-00 
056-50 C7Ee50 054.00 022200 
063-200 C8z2e00 060-¢00 036.00 
070200 CS3e00 066200 040200 
074.200 102.00 077.200 046200 
076250 1C€e50 085e00 0£1200 
079200 116200 089-00 055e00 
082250 127200 095.200 0€4-00 
9088.00 127.280 104.006 076200 
095200 14€.00 110e00 084-00 
106.200 154-00 11200 C8900 
121200 1€2200 113-50 0S&.00 
136-200 167200 114.50 1C1e00 
153200 168.30 117.250 106.00 
158e00 16€200 121-00 1120e00 
161-50 167230 127.00 116200 
166.200 166280 139-00 121-00 
171200 169 00 153.00 122.50 
183.200 168.00 165600 122200 
197.200 16€.00 179200 119200 
214.00 1€£206 193000 11560C0 
230-00 16€-00 213-00 104-00 
242.200 167290 229200 0922006 
250200 1€€.270 244-200 080e00 
253240 167220 257200 071200 
257200 169050 267200 O0€€e00 
265-00 169200 276200 062200 
280200 166200 289200 O£€e00 
291-00 1€2200 298200 049-200 
296200 159e30 307200 042.200 
298250 156290 318200 03£-00 
303200 159-00 333200 027200 
314200 193-00 347200 020¢00 
323200 148.00 360200 014200 
328200 142200 
331-00 1372e20 
335¢50 129220 
342250 121200 
349250 11€e00 
352250 113200 


III. COMPUTER PROGRAM DOCUMENTATION FOR PROGRAM-GRID 
1. Program Purpose. 


The purpose of this program is to determine the grid point array 
in the stretched curvilinear shelf coordinate system and appropriate 
scale factors needed for program SSURGE. The detailed grid is 
needed in order for the user to read off depths from an appropriate 
bathymetric chart of the shelf region at grid locations. 


2. Program Description. 


The program is written in FORTRAN IV language. This program 
interfaces between Programs CONFORM and SSURGE. It is assumed that 
the conformal mapping of the storm surge region has been completed 
to the user's satisfaction. The program GRID takes, in part, as 
input, the transformation coefficients and determines computing grid 
information of scale factors, grid point locations, and, at each 
grid point, the orientation of the €&-axis to the x-axis 


The program GRID is composed of 


MAIN Defines constants. Reads transformation coeffi- 
cients outputed from CONFORM. Reads water 
depths along a line near center of grid from 
the seaward boundary to the coast. Computes grid 
point locations, scale factors (uv, v , and FP), 
and, at each grid point, cos 6 and sine 6 
Writes computing grid information. 


SUBROUTINE XUT Writes information transferred into XUT. 


SUBROUTINE SHCOR Determines and writes grid point coordinates in 
X,y space and the distance in nautical miles 
between points. 


SUBROUTINE TRAN Computes the x and y coordinates of the trans- 
form-generated coastline and seaward boundary 
curve. 


SUBROUTINE TRAN1 Computes x(&,n) and y(é,n) 


SUBROUTINE TRAN2 Computes 9x/d— , dy/d— and 6 = tan-1 (22/28) 


SUBROUTINE CURV9 Contains the expansion curve Y = Z + B(x") 
, where A,B, and C are constants. The term 
Y is either Sp (units, nautical miles) or T 
(units, minutes). The term X is cither S* 
(units, nautical miles) or T* (mits, minutes). 
This subroutine computes Y and dY/dX given 
the coefficients and X . 
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SUBROUTINE CURV1 - Determines the necessary parameters to 
compute an interpolatory spline under 
tension through a sequence of functional 
values contained in arrays X2 and Y2. 


FUNCTION CURV2 - Interpolates the given curve, Y2 asa 
function of X2 , returning a value for 
Vi Salven gear 


SUBROUTINE CURV3 - Determines the necessary parameters to 
compute an interpolatory spline under 
tension through a sequence of functional 
values contained in arrays X2P and Y2P. 


FUNCTION CURV4 - Interpolates the given curve, Y2P asa 
function of X2P , returning a value for 
Va SAN CT Xane 


Type of Computer. 


The program GRID may be run on any computer with minimum 
core requirements of approximately 26K words (based on the 
present sample program appropriate to the Hurricane Carla 
surge simulation grid). GRID requires no auxiliary storage 
devices, peripheral devices, or magnetic tape input or output. 
No site-oriented computer plot routines are involved in the 
program. Approximately 25 minutes of machine time on a 
GE/635 is required for the sample program to determine the 
computing grid information. This time is based on the 
following pertinent program parameters: 


a) 150 transformation coefficients, BO or C, > 


b) 121 evenly spaced values of & for determining the 
transform-generated coastline arclength as a function 
OH EG 3 


c) 51 evenly spaced values of n for determining the arc- 
length along a particular isoline of &€ as a function 
Oi 8 


d) the computing grid of 47 € (or S*) lines and 15 n (or 
T*) lines; 


e) for determining the scale factor F , the area in x,y 
space of each quadrangle is approximated by using 4 
evenly spaced increments between €&€ isolines and 2 
evenly spaced increments between mn isolines. 
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For production runs, smaller sampling intervals might be required 
in) D, © eiatel Csjxscaailily, © - 


4. Input Data. 


Input data, other than constants defined in MAIN, are read in 
MAIN and are on IBM cards prepared according to the following list: 


(1) Card 1 


NMAX, NUMXI, Number of mapping coefficients, number of & 

NUMETA, DELSS, lines + 2, number of n lines, AS* in nautical 

SOPRA PDE ils miles; first value of S* in nautical miles, 

ND, NS AT* in minutes, number of depths, number of 
points in format 314, 3F5.1, 214. 


(2) Card 2 


BETA BZRO The value of 8 and By from the conformal 
mapping solution in format 2E14.7. 


(3) Card Group 3 


COB, COC - The NMAX cards containing the Fourier-type 
transformation coefficients, BL and C , 
in format 2E14.7. ” 


(4) Card Group 4 


SY Temporary storage for the ND values of the water 
depth (fathoms) along a line from the seaward 
boundary to the coast. This information is 
needed to evaluate the traveltime coordinate T. 


5. Mathematical Procedures and Program Limitations. 


4 

Information concerning the expanding grid procedure and the 
relations transforming &,n to S*,T* space is presented in 
Volume I of this report. The user is required to know the coeffi- 
cients of the expansion function 

S, = A + B(S*)” 

for each region of the curve where S. is arclength (nautical miles) 
along the transform-generated coastline. For the sample program, 
there are five regions of the expansion curve. Selecting AS* = 
6 nautical miles, the number of AS* intervals of each region and 
the value of dSp/3S* at the end points of each region, we can 
determine the coefficients of each region from three simultaneous 
equations derived from the constraints: 
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For region I, 176 nautical miles < S* < 236 nautical miles 
(10 intervals of AS* ) 


OSD 2 6205 Ae 8 


do* 


dSp 
35* 


S 
p 


Los) che | Sy 


S ale Se 
Pp 


For region IT, 236 nautical miles < 
(CMaintervallswore AS) 


OSD 3 '1..8 at S* 


dS* 


eS aye Se 


S aoe 


For region III, 260 nautical miles < 
(7 intervals of AS* ) ; 


3S* 


dSp 
3S* 


S 
wy 


Nos ge Ss 


1.0 ae §=6O 


302 nautical miles at S* 


= 176 nautical miles 


236 nautical miles 


236 nautical miles 


S* < 260 nautical miles 


236 nautical miles 


260 nautical miles 


260 nautical miles 


S* < 302 nautical miles 


= 260 nautical miles 


= 302 nautical miles 


= 302 nautical miles 


For region IV, 302 nautical miles < S* < 356 nautical miles 
(9 intervals of AS* ) 


SP 1.0 at S* 
SSP 1.0 aye }8=SR 
S = 302 nautical miles at S* 
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= 302 nautical miles 


= 356 nautical miles 


= 302 nautical miles 


For region V, 356 nautical miles < S* < 452 nautical miles 
(16 intervals of AS* ) 


SSP 1,0 at S* = 356 nautical miles 


asp 


aoe = Zoll at S* = 452 nautical miles 


S 


p 356 nautical miles at S* = 356 nautical miles 


A similar procedure is followed for the expansion function 
T=A+t BT*)° 


where T is the long wave traveltime (minutes) along the particular 
isoline of &. The program assumes that there are, at most, two 
regions of the expansion curve with the second region being AT = AT*. 
The expansion coefficients are determined by the program. If the 
user desires only one region (i.e., T = T* and AT* = total long wave 
traveltime-number of n grid intervals), the program computes AT*. 
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6. Flow Chart. 


MAIN - GRID 


Define constants 
NMAX, NUMXI, NUMETA, DELSS, DELTT, SSTRT 
ND, NS, XLAMDA, SLMAP, and XIDUM 


Read conformal mapping transformation 
coefficients, BETA, BZRO, COB(N), COC(N) 


Determine and write transform-generated coastline 

and seaward boundary curve coordinates for N4 evenly 

spaced ivalues ods: a (Om< ew) ae Store in ATXB , 
ATYB, and ATXA, ATYA. 


Determine and write coastline arclength 
‘Sp (nmi) as a function of &. 


Store Sp in X2 
Seo = a WA, 


CALL CURV1 
Fits interpolatory spline through Y2 
as a function of X2. 


Defines number of expansion regions (IRG) 

for Sp(S*) , coefficients A, B, and C and 

number of AS* grid intervals of each 
region (IIC). 
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Determine appropriate values of & for 
evenly spaced values of S* . 


For i = 1,2 -- NUMXI 


* 
S5 = SSSTRT + (i-1)*DELSS 


CALL CURV9 


Computes Sp for a given S* . 


FUNCTION CURV2 
Returns a value for & given Sp. 
Arrays X2 and Y2 cannot be al- 

, tered from the cast call of CURV1. 


Store THO YC1l 


E 


Determine and write coordinates and arclength 
Sn (nmi) along a particular isoline of € 
(=XIDUM) for NS evenly spaced values of n , 

=) <S i S Bo 
S 15 Ao 
n 
tim WZ. 


Store 


Store n 


CALL CURV 1 
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Read water depth (fathoms) at ND 
evenly spaced positions from the sea- 
ward boundary curve to the coastline. 


Determine and write the long wave traveltime 
T (min) and distance S, (nmi) along 
XIDUM for the ND positions. 


i ne 2P 
S fm | VWZIP 5 
n 


Store 


Store 


CALL CURV3 
Fits interpolatory spline through 
Y2P as a function of X2P. 


Define number of expansion regions (RG) for T(T*) , 

coefficients A, B, and C, and number of AT* grid 

intervals of each region (IIC). If IIC(1) = NUMETA-1, 
compute AT* 
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Determine appropriate values of jn 
for evenly spaced values of T*. 


For j - 1,2 .. NUMETA 


* 
T, = (j-1) * DELTT 


CALL CURV9 


Computes T for a given T* 


Function CURV4 


Returns a value for S, given 

T. Arrays X2P and Y2P can- 

not be altered from the last 
call of CURV3. 


Function CURV2 


Returns a value for nn _ given 

Sn AGkays exe andy 2yican 

not be altered from the last 
call of CURVI1. 


Store n in ETA. 


The values of & and n are now known. 
Determine the computing grid data. 
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CALL SHRCOR 
Determines and writes the grid point coordinates 
in x,y space and the distance (nmi) between 

grid points. 


Store x coordinate in X and 
Store y coordinate in Y for 
i = 2,3 .. NUMXI-1 
j = 1,2 .. NUMETA. 


Determine and write scale factor u 
X,y space-nmi). 


Store u in SX for i=1,2 .. NUMXI-2. 


(units 


Determine and write scale factor vw (units of 
X,y Space-time minutes). 
Store v in SY for j=1,2 .. NUMETA. 


Determine and write dimensionless scale factor F. 
Compute area in x,y plane of each quadrangle 
subdivided into IQUAD intervals between E 

lines and JQUAD intervals between n lines. 


Store area in X for i=1,2 .. NUMXI-1° 
1,2 .. NUMETA. 
=1,2 .. NUMXI-2. 
1,2 .. NUMETA. 


Store F anh SC Geo he 
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Determine and write cos 6 and sin 6 where 
@ is the orientation of the & axis to 
the x axis. 


Store 6 in Z for i=1,2 .- NUMXI 
j=1,2 .. NUMETA. 


Store cos 6 in X and 
Si -@ min SC stone ils .. NUMXI-2 
is 


1= 
j= NUMETA. 
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7. Description of Terms. 


Arrays 


All arrays except COB, COC, XI and ETA are reused throughout the 


program. 


The user is cautioned to consult each major program divi- 


sion for assessment of array contents. 


XO Ge Dimensioned NUMXI x NUMETA. 
X2,Y2 
X2P,Y2P 

ATXA,ATYA 

ATXB,ATYB Dimensioned the larger of N4, NS, ND, or IQUAD. 

Ne 135 (G ee 

TIC Dimensioned IRG (the number of Sp(S*) expansion 
regions). 

SX Dimensioned NUMXI-2. 

SY Dimensioned NUMETA. 

COB Dimensioned NMAX. Contains the conformal mapping 
transformation coefficients, Be 

Coc Dimensioned NMAX. Contains the conformal mapping 
transformation coefficients, ore 

XI Dimensioned NUMXI. Contains the values of & for 
determining the computing grid data. 

ETA Dimensioned NUMETA. Contains the values of n for 
determining the computing grid data. 

Constants 

NMAX Number of conformal mapping transformation coefficients, 
Bo @e G@ . 
n n 

NUMXI Number of computing grid € lines. With respect to the 
computing grid in Program SSURGE, there is an extra & 
line at each lateral end. This requirement results from 
the averaging procedure used in determining the grid 
data. 

NUMETA Number of computing grid n lines. This is the same 


number as in Program SSURGE. 
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DELSS 


DELTT 


SSSTRT 


ND 


NS 


IRG 


The value in nautical miles of AS* . This corresponds 
to DXI in Program SSURGE. 


Since in SSURGE the product, wAS* , is always computed in 
the surge equations, we determine DXI in meters such that 
u values are dimensionless, i.e., 


DXI = 1852 x AS * i cami | SLMAP —_2™ yy unit 
X,y unit 
The value in minutes at AT* . This corresponds to DETA 


in Program SSURGE. 


Since in SSURGE, the product, vAT* , is always computed in 
the surge equations, we determine DETA in meters such 
that v values are dimensionless, i.e., 


DETA = 1852 —4 
nm 


x AT* min. x SLMAP amt 


x,y unit 


The first value of S* in nautical miles. 


The number of water depths (fathoms) inputed from the sea- 
ward boundary to the coast for determining distance as a 
function of long wave traveltime. 


The number of points along XIDUM for determining n as 
a function of arclength Sh 3 


Horizontal extent of the mapped region in units of x,y 
space. 


The particular value of €&€ used in determining n(S_) 


The chart scale relating distance in nautical miles to 
distance in x,y units (i.e., nmi is equivalent to 5l 
units of length in x,y space). 


Acceleration due to gravity (feet-s~*). 


The number of points used in determining &€ as a function 
of arclength Sp. 


The number of expansion regions of SP(S*) or t(T*) 
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IQUAD The area in x,y space of each grid quadrangle is sub- 
and divided into IQUAD intervals between €& lines and 
JQUAD JQUAD intervals between n_ lines. 


8. Input and Output. 


The first card input to Program GRID gives the number of conformal 
transformation coefficients NMAX , the number of & lines NUMXI, 
the number of n Jlines NUMETA, AS* in nautical miles DELSS, the 
first value of S* in nautical miles SSSTRT, AT* in minutes DELTT, 
the number of water depths ND , and the number of points US _ used to 
establish n = n(S_) 


The second card gives 8 and B, and the next 150 cards give the 
mapping coefficients B, and C, . These 151 cards are the punched 
output from Program CONFORM. 


The remaining 31 cards are the depths picked off a bathymetric 
chart of the northwestern Gulf of Mexico. These depths are on a 
constant €&-line selected by the user. 


Expansion coefficients, provided by the user, appear as statements 
within the program after format 135. 


The reader must refer to Section III of Volume I of this report 
for an explanation of the parameters associated with the stretched 
shelf coordinate system for the Hurricane Carla surge simulation. 


Output from GRID consists of the transform-generated coastline and 
seaward boundary coordinates, the transform-generated arclengths along 
the coastline and seaward boundary, and for each of the five sections, 
dSp/aS* , S* , Sp , and €&€ are listed. Additionally, the 
transform-generated arclengths at values of n for evenly spaced 
increments of T* and at n values for constant increments of T 
along the chosen —-line, the n values along the specified €-line 
such that AT* is constant, and the traveltime and depths along the 
constant é-line are printed. 


The following output is required by Program SSURGE for simulating 
the Hurricane Carla surge. This includes the x,y coordinates of the 
grid intersections, the scale factors u and vw related to the trans- 
HORNACNOM Of & to Se ginvcl’ wy) wo We | megaeeemielhy, ie Seale iie- 
tor F , and the sines and cosines of theta giving the orientation of 
the €-axis to the x-axis. The program listing indicates when these 
may be punched or written on tape or disk for convenient input to 
SSURGE . 
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(s 

Cc 

Cc INPUT OATA 
150 47 15 €eC17600 beO 

Cc 


2ellO9S529EE01 $e5078648EE01 
— 403493873EE01-608176823EE01 
-301662011EE01 105427266EE01 
10242251 3EE0 1-101728220E600 
— 50 06365C2EEC0-10447071 8E-02 
106711 B64EE0 C= 106921 €334FE00 
— 305976 884EE 00-6073 30E€28E-01 
1041928E€0EE00 103667185EE00 
-101794944EE00-200S13£91E-01 
4e1190876E-01 1033577S3E-01 
—307759354E-01 48326 202E-03 
7eS716016E-02 £03212199E-02 
-1200800307E-01 106050344E-01 
—102900039E-C2=1 oS5SEB75E-01 
—5eSE35S7EE-C2 164022560E-01 
- 103795311 E-02-€ 03872108E-02 
-305084819F-02 926181509E-02 
606049823E-03-4 01274 305E-02 
—1LeS301L171E-C2 201339S555E-02 
Ge 3927 E01E-0 3-20 286E277E-02 
-101320465E-02 263924600E-02 
-100776871E-03=262€90670E-03 
-102749315E-CZ €032075S5E-03 
400667667E-03 104115754E-03 
—106674068E-03 40€354670E-04 
3099555 78E-0 4-206 3€ SEB TE-04 
-1288E€218SE-C3 5e682€825E-04 
900772255E-04 70€038S04E-04 
10559158 3E-0 3-4 09427928E=04 
Gel903310E-C4 30343SCS4E-04 
=101438246F-0 3-1 60257530E-04 
30e7952217E-04 167815754E-04 
-607125224E-04 9.487E230E-05 
4037 12E674E-04-507567649E-06 
-Se8352281E-C4 10470€E47E-04 
108351595E-04-3021 298S4E=-05 
—205129E44E-04 300412098E=-04 
= 100435642E= 05-2 04449992E-04 
100222874E-05 20632669 26E-04 
—5e8439003E-0€-10570C8E2E-04 
—102788127E-CE 105418€1 6E-04 
= 206134 267E-0 £=9 03255390E=05 
-1204402S52E-CE 100334477E=04 
= 500798231 E-6 €=502775377E-05 
-203834472E-0F 3024€6€107E-05 
20 6851102E-C £-206723854E-06 
~ 40311551 9E-0£-2003€ 34£2E-06 
206351927E-O5 1014225£3E-05 
— 203699234E-05 263080753E-06 
10060475 2E=-C £-2081665S6E-06 
-121909960E-0£F £0396716€2E-06 
7e96S51199E-06-703427355E-07 
—904350325E-0€ 26056501 4E=-06 
205921 445E=-06€-208028600E-06 
~ 30163786SE-CE £020276€6E-06 
800697 328E-07-3020493£3E-06 
~3elS5S61S70E-0E€ 20976EE27E-06 
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Si 


200821 353E-CE€-3010024£2E-07 
—20e7802158E-C6 7037£3£ 32E-07 
201309951 7E-CE-7e16EE7C2E-07 
—601223962E-C7 124500429E-06 
— 108451 7125-C7-10e11971C7E-06 
~303654205E-07 1eldDLES6SE-06 
Oe S384S1 SE~0E-405424002E-07 
—~3001203E€4E-C7 502354470E-07 
~3e6734064E-056—3e5E2570SE-07 
—1041408SCE-C7 4024405E5E-07 
703883 EB8SE-CE-10e2054ES8E-07 
—3022624E€6E-C7? 5e8860751E-08 
20 7314736E-C7? 1200232€4E-07 
—20802569SE-C7-202327544E-08 
2e23410S7E-07 126C852806E-08 
—904987154E-Ce& &.2818260E-08 
109075643F-0 E-€ 022565 73E-08 
—3e8SSIO06E-C2 706237684E-08 
102275 724E-CE-407253453E-CB 
—202807ESSE-OEF 50952221 76-08 
=407725860E-05-3eS1 76467E-08 
~1e 305SS312E-CE 42.9413841F-08 
1063466543E=-0 €-10e21 40226F-08 
— 301 784570E-CE-4.01840C8E-11 
2e5290277E-CE 102313703E-08 
20 8336831 E-C e-£ 054265 78E-09 
le 39OLSES4E-CE 402592CS21E-09 
“100413 04SE-CE 206€28°344F-095 
2094255S57E-05-329004823E-09 
~204943488E-CS 7205014S6E-09 
~201406560E-05-6 00237420E-09 
1e85C877SE-CS Ee3SS4S90E-09 
—307924576E-0 $-6092481S54E-09 
1eS02798S5SE-CS Ee0SEE2ZCEE-09 
—9e6309918E-10-205732490F-09 
-105450222E-0S 12536S5158E-09 
1e7949E78E-CS 3¢730S565E-10 
—204545E€29E-CS=-S eS SO4IELE-10 
109894453E-CS $e3920S30E-10 
—20e0184501E-09-€ 03272S67F-10 
103540433E-0S 70495€038F-10 
~10e18976S8E-0S—-2e61004€ 3E-10 
5e6070780E-10 GSe228SSHEE-11 
20872681 8E-1C 30045&C69E-10 
~Se8993701E-11-2e76E4ESOE-10 
9eS877297E-11 £olB3SE14E-10 
-109201876E-10-308510567E-10 
6e7366526E-11 364544254F-10 
~8eS2129S1E-11-201782821E-10 
—306849378E-12 12826C0446E-10 
~30e18S53177E-12-00 295704 3E-11 
-5e8541698E-11 6GelfESB99F=-11 
Se4400211F-11-720616809E-12 
—7064321740E-11 203844000E-12 
5Se4853889E-11 1e7E€02C39E-11 
—5e9319842E-1] 1-101 7E3087E-12 
3eSO82522F-11 7e7SEESG7T7TE-12 
—303595148E-11 304384777E-12 
109693600E-11 368907201E-13 
~2006267S8E-11 SeE511549E-12 
100770427E-1 1-1 0€112022E-12 
—100974893E-11 €e770SE36E-12 
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30 7490653E-12-500353240F-12 
-306€CS5S470E-12 801811602E-12 
= 405852351E-1 3-£ 0 946SE72E-12 
= 704799452E-12 7643025436-12 
= 10297923€E-12-SelE12560E-12 

90 3319ES4E-14 503514749E-12 
-8e512744B8F—1 3-3 0465ES60E-12 
—4e6265112F-12 202230€44E-12 

2074291 5AF-1 2-10 365E647F=12 
—105414227E-12 Ee 782235HE-13 

1Lef202911E-12 2053572546E-13 
— 1054622 74E-12-2 08ES23E4E 13 

Le 42H7FF3E-12 60296€02464F-13 
= 104235914F-1 2-1 69242733E-13 

8e0951983E-12 204162215 =-13 
—6074564SEE-12 161457465E-13 

206952090E-13-165477038E-13 
—10687.7584E-12 460572501E-13 
= 909264500F=1 4-3 0736£145E=13 
609077629E=14 “40 7001E€51E=13 
- 102943830E-12-2040€C812E-13 

10 7026387E-14 20511€468E-13 

102898776E-15=102325175E-13 
-100197686E-13 $e0332572E-14 

90003627 4E- 14-6 03540502E-15 
= 10253069SE-12-1 00823757E=15 

100618124F-12 4643203E8E-14 
—101324192E-12-163624101E-14 

7e0773E€32E-14 20450EE13E-14 
—602031906E-14 30025€047E-15 

304 432727E-1 4-1 01037802E-14 
-1097S2800E-14 3e7E0CE72E-14 
956 
636 
53e¢ 
460 
390 
32e 
300 


IV. COMPUTER PROGRAM DOCUMENTATION FOR PROGRAM-SSURGE. . 
1. Program Purpose. 


The purpose of the program is to numerically simulate the storm 
surge in orthogonal curvilinear coordinates with the vertically 
integrated form of the quasi-linear long-wave equations. 


2. Program Description. 


The program is written in FORTRAN IV language. It is assumed 
that the conformal mapping of the region under investigation has 
been completed. The transformation coefficients for three regions 
of the continental shelf of the Gulf of Mexico and two regions of 
the eastern seaboard are provided in Appendixes A and E in Volume I 
of this report. 


An interfacing program is required which inputs the coefficients 
and generates a curvilinear computing grid to the user's satisfac- 
tion. The output from the program (and, in part, also the input to 
SSURGE) must be the scale factors, grid point locations, and, at each 
grid point, the orientation of the €-axis to the x-axis. 


The Program SSURGE is composed of: 


MAIN Defines constants and interfaces the 
subroutines. 


SUBROUTINE ZERO Initializes all arrays to zero. 
SUBROUTINE FIELD Reads data and writes the water depth field 


relative to mean sea level, the wind 
field parameters and the storm positions. 


SUBROUTINE WINDF Calculates the wind and atmospheric pressure 
fields. 


SUBROUTINE ELEV Computes the water tevel anomaly, H . 
SUBROUTINE FLUX Computes transports, Qox and Qn ¢ 
SUBROUTINE DRAW1 Outputs H and vertically averaged water 


velocities, Qou/D and Qx/D , 


at hourly time intervals and saves the water 
level anomaly at prescribed grid locations 
for output at program completion. 
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SUBROUTINE METER Calculates and saves the vertically averaged 
water velocities at prescribed grid locations 
for output at program completion. 


SUBROUTINE HUV Outputs the saved information of the simu- 
lated hydrographs, simulated current meters, 
and observed water levels at program completion. 


Type of Computer. 


The program SSURGE may be run on any computer with minimum 
core requirements of approximately 30K words of memory (based on 
the present sample program appropriate to the Hurricane Carla 
surge simulation). The program has been executed successfully on 
the IBM 360, CDC/6600 and 7600, and GE/635. The present sample 
program requires no auxiliary storage devices, peripheral devices 
or magnetic tape input or output. No site-dependent computer 
plot routines are involved in the program. Approximately 14.4 
minutes of machine time on a GE/635 is required for the sample 
program to complete 66 hours of surge simulation. 


Input Data. 


Input data, other than constants defined in MAIN, are read 
in SUBROUTINE FIELD. These data are on cards prepared according 
to the following list. 


(1) Card 1 


NT1 - Number of cards (max.50) containing on each TIM, ROT, 
RAD, VRMAX and PZRO (format I5). 


(2) Card Group 2 


NT1 cards with each card containing values of 


a) TIM The time in hours at which the hurricane wind 
and atmospheric pressure parameters are 
recorded (format F10.1). 


b) ROT The angle in degrees between the direction the 
storm is moving and the region of maximum winds 
(format F10.1). 


c) RAD The distance in nautical miles from the storm 


center to the region of maximum winds 
(format F10.1). 
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(3) 


(4) 


(5) 


(6) 


(7) 


d) VRMAX Maximum observed windspeed in knots (format 


F10.1). 


e) PZRO Atmospheric pressure in millibars of the storm 


center (format F10.1). 


Card 3 


NT2 - Number of cards (max.150) containing on each 


TIMPOS, XPOS and YPOS (format I5). 


Card Group 4 


NT2 cards with each card containing values of 


a) TIMPOS The time in hours at which the hurricane 


position is recorded (format F10. 


A) 


b) XPOS The x-coordinate in units of x,y space of 
the hurricane center (format F10.1). 


c) YPOS The y-coordinate in units of x,y space of 
the hurricane center (format F10.1). 


Card Group 5 


GRID2 The fluid depth in fathoms relative to 
level along each column, i=1,2..IM, is 
a nested do-loop for j=1,2..JM (format 
The program will zero those values for 
prior to computations. The depth data 


mean sea 
read with 
11F7.2). 
even itj 
are posi- 


tive numbers which the program converts to nega- 
tive values (in meters) to be consistent with 


the coordinate system. 


Card Group 6 


S The values of the dimensionless scale factor relating 


the (x,y) plane to the (&€,n) plane are read along each 
column, i=1,2..IM, with a nested do-loop for j=1,2..JM 


(format 5E14.7). 


Card Group 7 


DSDXI The values of the dimensionless scale factor, u, 


transforming &€ to S* are read with 
per card for i=1,2..IM (format E14.7). 


SiG 


one value 


(8) Card Group 8 


DTDET Values of the dimensionless scale factor wv 
transforming n to T* are read with one 
value per card for j-1,2..JM (format 2X,E14.7). 


(9) Card Group 9 


HOBS1 The values of the observed hourly water level 


in feet at grid location (IH1, JH1) are read 
with 19 values per card in format F4.1. 


(10) Card Groups 10 through 14 


HOBXk The values of the observed hourly water level 
in feet at grid location (IHk, JHk) are read 
sequentially with the same format as above. 


(11) Card Group 15 


YOSo 2 LAL The paired x,y coordinates (units in x,y 
space) of the computational grid points are 
read along each column, 1-1,2..IM, with a 
nested do-loop for j=1,2..JM (format 10F7.2). 


(12) Card Group 16 


COSG, SING The paired values of the cosine 6 and sine 6 
where 6 is the angle between the € and x 
axis are read along each column, i=1,2..IM, with 
a nested do-loop for j=1,2..JM (format 10F8.5). 


A computer printout of the sample program and data cards are given 
later in this section. 


5. Mathematical Procedures and Program Limitations. 


The storm surge equations, solutions, and algorithms are described 
in the text of Volume I of this report. 


Basically, the model utilizes a centered difference, leapfrog 
analog of the vertically integrated, quasi-linear form 
of the long-wave equations. The algorithm treats the time 
dependency explicitly and employs a computing lattice in which 
the transports, Qox and Qa >» are computed at the same location 


but are staggered in time and space with respect to the water 
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level anomaly. The program assumes that seabed scouring does 
not occur. No program error messages or consistency checks are 


incorporated in this version. 


Flow Chart. 
MAIN - SSURGE 


Define constants 


CALL ZERO 
Initializes all arrays to zero 


CALL FIELD 
Reads data and writes depth field, 
wind field parameters and storm positions 


CALL HUV 
Output saved information of 
simulated hydrographs, simulated 
current meters and observed 
water levels. 
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CALL WINDF 
computes wind and atmospheric pressure 
fields WIND(i,j) , P(i,j) 


CALL ELEV 
computes water level anomaly GRID1 (i,j) 


CALL DRAW1 _ 
outputs H, Qcox%/D and Qrx/D at hourly time 
intervals and saves information for output 
at program completion 


CALL FLUX 
computes transports, GRIDI. (i,j) and GRID2(i,j) 
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Glossary and Description of Terms. 


Arrays 


GRID1 Dimensioned IM x JM. Contains H values in meters 


at oad a) and Qo in meters squared per second 
at even itj 


GRID2 Dimensioned IM x JM. Contains fluid depth values, 
Dy » in meters at odd itj and Qps in meters squared 
per second at even itj 


S Dimensioned IM x JM. Contains dimensionless scale 
factor, F, determined from the conformal mapping of 
(x,y) space to (&,n) space. 


DSDXI Dimensioned IM. Contains the dimensionless scale 
factor, u, transforming §& to S* 


DTDET Dimensioned JM. Contains dimensionless scale 
factor, v , transforming n to T* 


XX Dimensioned IM x JM. Contains the x-coordinate in 
units of (x,y) space of the grid point locations. 


YY Dimensioned IM x JM. Contains the y-coordinate in 
units of (x,y) space of the grid point locations. 


COSG Dimensioned IM x JM. Contains cosine values of 0 
relating the orientation of the &-axis to the 
X-axis at the computational grid points. 


SING Dimensioned IM x JM. Contains sine value of 6 
relating the orientation of the &-axis to the 
X-axis at the computational grid points. 


WIND Dimensioned IM + 1 x JM. Contains values of the 
windspeed in meters per second. At a transport 
computational grid point (i,j), the S*-component | 
is stored in WIND at (i,j) and the T*-component 
is stored at (itl,j). 


P Dimensioned IM x JM. Contains values of Hg (the 
hydrostatic elevation in meters of the sea surface 
corresponding to the atmospheric pressure anomaly) 
and stored at odd itj. 


HOBSk Each array is dimensioned 67. Contains the observed 


en.2. 6 hourly water levels in meters at grid locations 
a (IHk, JHk). 


4I 


ROT 


TIMPOS 


XPOS 


YPOS 


Constants 


IM 


JM 


NMAX 


Each array is dimensioned 300. Contains values 
of the computed water level H in meters at 
grid locations (IHk, JHk). 


Each array is dimensioned 300. Contains values of 
the vertically averaged water velocity, Qc»/D , in 
meters per second at grid location (ITk, JTk). 


Each array is dimensioned 300. Contains values _ 
of the vertically averaged water velocity, Qyx/D , 
in meters per second at grid location (ITk, JTk). 


Dimensioned 50. Contains the time in hours at 
which the wind and atmospheric pressure field 
parameters (ROT, RAD, VRMAX and PZRO) are 
recorded. 


Dimensioned 50. Contains the angle in degrees 
between the direction the storm is moving and 
the region of maximum winds. 


Dimensioned 50. Contains the distance in nautical 
miles from the storm center to the region of 


- maximum winds. 


Dimensioned 50. Contains the maximum observed 
windspeed in knots. 


Dimensioned 50. Contains the atmospheric 
pressure in millibars of the storm center. 


Dimensioned 150. Contains the time in hours at 
which the storm position is recorded. 


Dimensioned 150. Contains the x-coordinate in 
units of (x,y) space of the storm center. 


Dimensioned 150. Contains the y-coordinate in 
units of (x,y) space of the storm center. 


Number of grid points in the S* direction. 
Number of grid points in the T* direction. 


Maximum number of time steps. 
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INC 


DXI 


DETA 


DELT 


F 
THk , JHk 
keig25 56 


ITk, JTk 
k=1,2..6 


CORIO 
PHI 


PINF 


Comments 


Number of time steps between saving of surge 
results in HYDk, UCMk and VCMk. INC must be 
an even integer number. 


Grid increment in meters in the S* direction. 


Grid increment in units of meters in the T* direc- 
tion. The units of DTDET-DETA (that is, v-AT*) are 
in meters. 


Time increment in seconds. DELT must be less 
than that required for numerical stability and, 
also, be an even integer multiple of 60 minutes. 


Acceleration due to gravity (=9.8 meters per 
second squared). 


Dimensionless seabed drag coefficient (=0.0025). 


Grid point location for saving the computed water 
level. The sum of the indexes must be odd. Index 
THk cannot exceed IM. Index JHk cannot exceed JM. 


Grid point location for saving the vertically 
averaged water velocities. The sum of the 
indexes must be even. Index ITk cannot exceed 
IM. Index JTk cannot exceed JM-1. 


Coriolis parameter (=6.70875 x 10 © second 1) 
for latitude: 27s 256252! Ne 


Wind ingress angle in units of degrees. 


Far field atmospheric pressure in millibars 
(=1016 millibars). 


1. A symmetric analytical hurricane wind field representation 
as given by C. Jelesnianski (1965, A numerical calculation 
of storm tides induced by a tropical storm impinging on a 
continental shelf, Mon. Wea. Rev., 94, 379-394) is 
employed in the surge model). 


2. The wind stress coefficients, K, and Ky , are 1.1 x 107° 
and 2.5 x 107-© and are defined in SUBROUTINE FLUX. 


3. Constants YRANGE, THIT, XHIT and YHIT are not used in 
this program version. 
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8. Input and Output. 


Data statements in Program SSURGE establish the number of computa- 
tional points, number of time steps, output interval, values of the 
grid and time steps, acceleration due to gravity, bottom stress coeffi- 
cient, and locations of the simulated hydrographs and current meters, 
and the Coriolis parameter corresponding to latitude 27° 23.232' N. 


Card input provides the number (NT1) of hurricane description sets, 
followed by NT1 cards giving the time (in hours after start of com- 
putations) and the three required storm parameters. These are suc- 
ceeded by one card giving the number (NT2) of hurricane positions to 
be used, followed by NT2 cards providing the time and storm center 
positions in the original arbitrary Cartesian grid (see Volume I of 
this report). NT1l need not be equal to NT2, nor must the observed 
storm parameters and positions coincide in time. 


For the Carla surge computations, the hurricane characteristics are 
stipulated at 6-hour intervals for the first 18 hours and at 3-hour 
intervals for the remainder of the 66-hour prototype time simulation. 
Note the radius to maximum winds are in nauttcal miles, the maximum 
winds are in knots, and the central pressures are in mtllibars. Due 
to the erratic movement of Hurricane Carla, the coordinates of the 
center of the storm are specified at hourly intervals, except for a 
single 6-hour interval spanning the end of the simulation. These 
coordinates are specified in x,y space (100 units = 219 kilometers). 


The depth field, in fathoms, is introduced followed by the array of 
scale factors, S . The array of scale factors, DSDXI (an alias for 
u ) is read along the specified row and the scale factor array DTDET 
(an alias for wv) is specified along the chosen column. 


The following six arrays are the observed hourly water levels (in 
feet) from Padre Island (HOBS1), Aransas Pass (HOBS2), Port O'Connor 
(HOBS 3), Pleasure Pier, Galveston (HOBS4), Mud Bayou (HOBSS), and 
Sabine Pass (HOBS6). These data are not necessary for any phase of the 
surge calculations and may be omitted. 


The last two arrays input to SSURGE give the coordinates (XX,YY) 
of the grid points in units of x,y space, and the sines and cosines 
of theta at the computational points. 


All input not in the MKS system of units is converted internally 
to the MKS system. 


The storm parameters, hurricane center coordinates, and depth field 
are printed out immediately following input. At hourly time intervals, 
the water level anomalies and the depth averaged velocity component 
fields are printed. The six simulated hydrographs and current meters 
are printed out with their positions and the corresponding observed 
water levels at the completion of the surge simulation. 
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